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We present microscopic wave functions of spin singlet Mott insulating states and nematic 
Mott insulating states. We also investigate quantum phase transitions between the spin sin- 
glet Mott phase and the nematic Mott phase in both large-iV limit and small-iV limit (N 
being the number of particles per site) in high dimensional bipartite lattices. In the mean field 
approximation employed in this article we find that phase transitions are generally weakly first order. 
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I. INTRODUCTION 

The recent observation of correlated states of bosonic 
atoms in optical lattices has generated much interest 
As known for a while, when bosons in lattices interact 
with each other repulsively, they can be localized and 
form a Mott insulating state instead of a condensatei^ 
This phenomenon has been observed in the optical lattice 
experiment. By varying laser intensities of optical lat- 
tices, Greiner et al. have successfully investigated Mott 
states of spinless bosons by probing spin polarized cold 
atoms in optical lattices with a large potential depthi^ 

We are interested in spin correlated Mott insulating 
states of spin-one bosons, especially spin-one bosons with 
antifcrromagnetic interactions. Some aspects of spin cor- 
related Mott insulating states were investigated recently. 
For an even number of particles per site, both spin sin- 
glet Mott insulators and nematic Mott insulators were 
found in certain parameter regimes, while for high di- 
mensional lattices with an odd number of particles per 
site only nematic insulating states were proposed^ In 
one-dimensional lattices, it was demonstrated that for an 
odd number of particles per site, Mott states should be 
dimerized valence-bond-crystals, which support interest- 
ing fractionalized quasi-excitations£ Effects of spin corre- 
lations on Mott insulator-superfluid transitions have been 
studied and remain to be fully understood^ 

In this article, we analyze the microscopic structures of 
spin singlet Mott insulating states (SSMI) and nematic 
Mott insulating states (NMI). We study, quantitatively, 
quantum phase transitions between these two phases in 
high-dimensional bipartite lattices. In the mean field ap- 
proximation, we demonstrate that for an even number 
of particles per site, the transitions are weakly first or- 
der. We should emphasize that results obtained in this 
paper are only valid in high dimensions. In one dimen- 
sional lattices, nematic order does not survive long wave 
length quantum fluctuations; detailed discussions on low 
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dimensional Mott states for both even and odd numbers 
of particles per site are presented in Ref. 6. 

The organization is as follows. In section II, we present 
the general setting for the study of spin order-disorder 
quantum phase transitions. In section III, we present 
mean field results on the quantum phase transitions in 
both small N and large N limits. In section IV, we dis- 
cuss issues which are to be understood in the future. 



II. ALGEBRA AND SETTING 

A. The microscopic Hamiltonian in the dilute limit 

The microscopic lattice Hamiltonian we employ to 
study spin correlated states of spin-one bosons is: 

^microscopic — ^ (^fc.m^; m ^ h.C.) 

(kl) 

+ ^lrn^rnU P (k,l)i>l ml i> ltml (1) 
k,l 

k,l 

Here tp k is the creation operator of a spin-one particle 
at site k with spin-index m = 0,±1. (kl) indicates that 
the sum should be taken over nearest neighbors and 5 7 
(7 = x, y, z) are spin-one matrix operators given as: 

1 f° 1 °\ 1 f° ~ l °\ ( l °\ 

S x = - r \ 10 1 S v = - r \ i -i] S z = . 

v 2 \0 1 0/ v 2 \0 i 0/ \0 0-1/ 

U p (k,l) and U s (k,l) are, respectively, spin-independent 
and spin-dependent interaction parameters between two 
bosons at site k and I. 

In the dilute limit, which is defined as a limit where 
pa? <C 1 (a is the scattering length and p the average den- 
sity), atoms scatter in s-wave channels. For two spin-one 
atoms, the scattering takes place in the total spin S = 0, 2 
channels, with scattering lengths ao,2- Interactions be- 
tween atoms can be approximated as spin-dependent con- 
tact interactions^ In the lattice model introduced here, 
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calculations yield 

U»(k,l) = E c 5 kl and U s (k, I) = E s 5 kl . (2) 
The parameters E c and E s are given by: 

A-kK 2 p{2a 2 + a ) . „ 47rn 2 p(a 2 - a ) . 

& R = „, C, ii„ = rTTTTT C, (3) 
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where iV is the average number of atoms per site, M is 
the mass of atoms and c is a constant. 



Of particular interest is the singlet creation operator 
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g*W*,a " ^g(<0<0 ~ 2 ^,l^,-l)- ( 10 ) 

We find the following properties for this operator: 

- [S^,^,a] = 0; (Ha) 

[ipk^ipk^^lp] = 2Skiipk,p; (Hb) 

[tpk,a^k,a,^j^tplp} = 5fcj(4p fe + 6). (11c) 



t ./,t 



B. Algebras 

For the study of spin correlated states in lattices, it is 
rather convenient to introduce the following operators: 



= ^=w»I,-i-vii) 



v> 



(4a) 
(4b) 

= 4,0 ( 4c ) 
where fc again labels a lattice site. In this representation: 



t - J_ 

k, v ^ 

t 



(5) 



a, /?, 7 € {.t, y, 2}. The density operator can be expressed 
in a usual way: 

Pfe = 4>l n i>k,i- (6) 

Consequently the Hamiltonian is given as: 

Wlatt=-*X (Via V>i,a+ n.C.) + ^ c p| + E s S|-^ (7) 
(fci) ft 

where we have introduced the chemical potential p; S 2 . 
is the total spin operator Sk. a Sk,a- 

^k,a (fit = x, y, z) are bosonic operators obeying the 
following commutation relations: 



[i>k,a,1pl,p] = bl>k, a iVl,p\ = °> bPk^^lf}) = falScc/3- (8) 

Taking into account Eqs. 5,6,8, one can verify the fol- 
lowing algebras: 



[S£,sf] 


= 5 kl ie a ^ k „ 
= 6 kl ie a ^l n 
= d kl ie a ^Sl 


(9a) 
(9b) 
(9c) 


[Pfc,^i,a] 


= —8kl1pk,a 

= <wi, Q 


(9d) 
(9e) 




= 


(9f) 



C. The on-site dynamics 

The total spin operator can be expressed as: 

si = Pkifik + 1) - i>l, a Tpt a i>k,(3iPk,0- (12) 

So, eigenstates of the total spin operator have to 
be eigenstates of the "singlet counting operator" 

V'LViaVMV'M- 

Defining the state \&£ such that: 

we find that wave functions of these eigenstates are: 

where C is a normalization constant. From Eq. Illcl it 
follows that: 



4,A,«Mwn, m = (4m(n - m) + 2m)n,m (15) 
Using that pk^k,m = n ^k,m we derive: 

Stn.m = (n - 2m) (n - 2m + l)¥j? ro . (16) 



So 5fc = n — 2to. Now if n is even, S k is also even and 
when n is odd, S k is odd too. For an even number of par- 
ticles per site N the states labeled by S k = 0, 2, 4, . . . , N 
are present, whereas for an odd number of particles per 
site S k = 1,3,5, ... ,N are allowed. This reflects the ba- 
sic property of the many body wave function of spin-one 
bosons, which has to be symmetric under the interchange 
of two particles. 

Solutions for spin correlated condensates with finite 
numbers of particles were previously obtained^; in the 
thermodynamical limit, these states evolve into polar 
condensates&i2*ii Also there, two-body scatterings were 
shown to lead to either " antiferromagnetic" or " ferromag- 
netic" spin correlations in condensates. Spin correlated 
condensates have been investigated in experiments^*^ 



D. The effective Hamiltonian for Mott states 

In the limit when t <C E c , atoms are localized and 
only virtual exchange processes are allowed. An effective 
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Hamiltonian in this limit can be derived in a second order 
perturbative calculation of the Hamiltonian in Eq0 



S 2 / 



(kl) 



(17) 



Here J e 



2^-. In deriving Eq. 1171 we have taken into 
account that E s <C E c . 

To facilitate discussions, we introduce the following op- 
erator: 



Qk,a/3 = Ipt a ^k,f3 ~ ■zSaplpt ^k tl , (18) 



whose expectation value 



Q 



(Qap) 



(19) 



is the nematic order parameter. The reference state 
-01-of = IIa; — ° / N ° — 10) is a maximally ordered state. 



Choosing n = e z , we obtain: 

(Qc/3)rcf = N 



-| 
-i 

o I 



(20) 



Q varies in a range of [— ~ , 1] . 

In terms of the operator Q a p, the effective Mott Hamil- 
tonian can be rewritten as (up to an energy shift): 

n^ = E a ^Tt[Q k Q k -Q k Q[]-J m ^T[[Q k Qt]. ( 21 ) 



(k-i) 



Finally we define 



n 



zJ c 
E„ 



(22) 



as a dimcnsionless parameter, which can be varied con- 
tinuously; z is the coordination number of lattice. 



III. PHASE TRANSITIONS BETWEEN SSMI'S 
AND NMI'S 

A. Two particles per site 

In the case of two particles per site , the on-site Hilbcrt 
space is six-dimensional, including a spin singlet state 



|S = 0,S z = 0) = ^|0>, 



and five spin S = 2-states 



(23) 



(24) 



where Q v ^ is a symmetric and traceless tensor with five 
independent elements. All states in the Hilbert space are 
symmetric under the interchange of bosons; as expected, 
the states \Q v t) arc orthogonal to IS = 0, S z = 0). It is 
convenient to choose the following representation of Qr^'. 



(25) 



with the director n as a unit vector living on S 2 . States 
defined by the director n form an over-complete set in 
the subspace spanned by five S = 2 states. 

When the hopping is zero, one notices that the Hamil- 
tonian in Ea ll7l commutes with S 2 ; the ground state wave 
function is 



I*) 



~[\S = 0,S z = Q) k . 



(26) 



On the other hand, when E s goes to zero, the Hamilto- 
nian commutes with Tr[Qfc,a/3<3z,/3a] and the ground state 
wave function can be confirmed as: 



*) = II \I\\Q^))u + ^=\s = o,s z = o) 



(27) 



E. The range of the physical parameters 

From Eq. [3] it is clear that E s and E c depend on the 
density, number of atoms, the mass of atoms and scat- 
tering lengths. However, their ratio depends only on the 
scattering lengths. According to current estimates^i^, 
for sodium atoms this ratio is given as jf- ~ 9 ■ 10~ 2 . In 
this paper, we are interested in the limit E s <C E c . 

The parameter t can be varied independently by chang- 
ing the depth of the optical lattice. A wide range is exper- 
imentally accessible; one can vary from the regime where 
i ^> E c to a regime where t <C E s . We limit ourselves to 
Mott states (i <C E c ), where all bosons are localized, but 
the ratio fj can have arbitrary values. 



for any choice of the director n. 

To study spin nematic or spin singlet Mott states at 
an arbitrary fj, we introduce a trial wave function which 
is a linear superposition of singlet states and symmetry 
breaking states: 



= n cos0 |S = 0,S, = 0)k + sme\Q vi (n)) k . (28) 



Here 9 is a variable to be determined by the variational 
method. 

A straightforward calculation leads to the following re- 
sults: 
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FIG. 1: Energy (measured in units of E s ) versus Q for various 
fj for N = 2. Curves from top to bottom are for fj = 0.97, 
0.99, 1.0, 1.02. 



E{9) = 



z 2 



= 6£ s sin 2 9 - zJ cx - ((2\/2cos 9 sin 6> + sin 2 9) 
3 



(29) 

2 



Q = (a/2 cos 6* sin ( 



sin 2 6> 



(30) 



In terms of Q, the energy can be expressed as: 



which for Q -C 1 can be expanded as: 



2 



(31) 



The cubic term leads to a first order phase transition 
in the mean field approximation, which is similar to the 
situation in classical nematic liquid crystals^ 

In FIG^ the Q-dependence of energy is plotted for 
various fj in the vicinity of a quantum critical point (mean 
field). For fj < 0.985, the energy has only one minimum 
at Q = and correspondingly the ground state is a spin 
singlet Mott state. When 0.985 < fj < 1.0, in addition 
to the global minimum at Q — 0, there appears a local 
minimum at Q > 0, which represents a spin nematic 
metastable state. When fj > 1.0 the solution with Q > 
becomes a global minimum and the solution at Q = is 
metastable; consequently the ground state is a nematic 
Mott state. For fj > |, the solution at Q — becomes 
unstable ; but an additional local minimum appears at 
Q < which we interpret as a new metastable state (not 
shown in FIG. [J). 

The evolution of ground states as fj is varied, is sum- 
marized in FIG. [3 As is clearly visible, the phase tran- 
sition is a weakly first order one. The jump in Q at the 
phase-transition (fj = 1.0) is equal to ^. 
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FIG. 2: (Color online) The nematic order parameter as a 
function of fj for N = 2. The phase transition takes place at 
fj = 1. Data along the black lines represent ground states; 
the red (light) lines are for metastable states. Spheres with 
double-headed arrows are introduced to represent ordering in 
director n defined in Eq. 25 in different Mott states. In spin 
singlet states, the director n is uncorrelated; in rod-like ne- 
matic states, the director n is ordered and in disk-like states, 
the axis of the easy plane of the director n is ordered. 



It is worth emphasizing that a positive Q corresponds 
to a rod-like nematic state; for Q = 1 the state is micro- 
scopically given by: 



V2 



|0>. 



(32) 



A solution with negative Q indicates a disk-like ne- 
matic state; the microscopic wave function is 



v4v4-™Mlo>. 



(33) 



at Q = — i. For n = e z the wave functions in 
Eq.32,33 become ^M|0) and ± (^M + |0) re- 

spectively. 

We have also tried a five-parameter variational ap- 
proach, taking into account the full on-site Hilbcrt space. 
In a slightly different representation we write the trial 
wave function as: 

|*> = Yl ( c xx\xx) k + c V y\yy)k + c zz \zz) k 
k 

+c xy \xy) k + c xz \xz) k + c yz \yz) k ) (34) 

Here \aa)k = -j= ^.q^.q |0) (no summation) and 

\a(3)k = ip k a ip\ a|0). This results in the following ex- 
pression for the energy: 
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E = E s [4(cL 

-2 Jox [6((4 



-■xx^yy 



ci, + 4) + 4(4 + 4 



i i/ 2 2 , 2 2 , 2 2 \ 
~*~^\ c xx c yy "r ' c yy c zz I 



^yy^zz) 



" 4 ) 



„2 „2 



"sir 

^l^( C 2;a; C 2;i/ "I" C a;a: C a;z C yy C xy ^ 

. 2 2 | 2 2 \ 
~r c zz c a:z c zz C !/zj 

~*~®\ c xy c xz i c xy c yz ' c iz c ^zJ 
_ F8V^(c ;c:r -|- Cyy -\- C zz )c xy C xz C y , 



2 2 

yy yz 



^/-* 2 _L ^ 2 ^ 2 

>^\ c xx c yz ' c yy c xz 



" 2 ) 

-'xyJ 



~\-^>{C XX CyyC xy + C XX C ZZ C XZ + CyyC ZZ C xz )\ 



(35) 



The conclusions are almost the same and summarized 
below: 

i) For fj < 0.985. the only minimum is at c xx = c yy = 
c Z z = ■jg, c af3 = for a ^ /3. 

ii) At fj = 0.985 additional local minima appear. 

iii) At fj = 1 a first order phase transition takes place. 

iv) For I > fj > 1, the global minimum is at Q > 0, but 

the Q = 0-solution remains to be a local minimum. 

v) At fj = § the solution at c xx = c yy = c zz = -4= 
becomes unstable. 

vi) However, the disk-like Q < 0-solution appears in this 
case as a saddle point. 



The Hamiltonian in Eq. 7 can be mapped to a Con- 
strained Quantum Rotor Model (CQR), describing the 
dynamics of two unit vectors (n, e lx ) on a two-sphere 
and a unit circle: 



W C qr = -t^n/fc-n/ cos(xk-Xi)+^2E s S 2 k +E c p 2 k -p k fi 

(kl) k 

(41) 

t = Nt. The CQR-modcl has been introduced to study 
spin-one bosons in a few previous works and we refer 
to those papers for detailed discussions Ji&lL For Mott 
states the effective Hamiltonian can be found as: 



H = E S J2SI- Jcx^K-n,) 2 ; J ex = — , (42) 



(kl) 



and we define 77 = ^ 



In general, we choose the on-site trial wave function to 



be: 



tp(n k ) = G G exp ^(n fc ■ n ) 2 



(43) 



B. Large N limit: An even number of particles per 

site 

For a large number of particles per site, it is convenient 
to introduce the following coherent state representation: 



N+SN 



|n,x) = 



V26N, 



^2 exp(-imx) 



Mi)' 



m=N-SN 



y/2(m-l)\ 



|0) (36) 



where the director n is again a unit vector on S 2 given 
by (cos cj) sin 6, sin (j> sin 9, cos 6*). In this representation 



P 
S 

s 2 

Qaf3 



. d 

£>xk 



= in x 



d_ 
dn 

1 d 



1 8 2 



sin 



N n„n 



la/3 



(37) 
(38) 
(39) 
(40) 



Co is a normalization constant. When a — *• this yields 
an isotropic state loo( n fc); which indicates a spin singlet 
state. When a — > +00, is localized on the two-sphere 
in the vicinity of n , representing a rod- like nematic state 
and when a — > —00, lies in a plane perpendicular 
to no, corresponding to a disk- like spin nematic state. 
Moreover this wave function has the following property: 
ip{— rife) = i/>(nfe), as is required for an even number of 
particles per site^ 

Choosing no = e z this gives: 



ip{4>k,0k) = C CT exp 



cos 9k 



(44) 



The expectation value of the Hamiltonian in this state 
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Ea = Es Ml + 3eVH 

I 4 2 2ViErfi^| 



12e 2<T cr - 4e< T 0V Mt 3 + 2cr)Erfi v /|cr| + tt(3 + 4(cr + cr 2 ) )ErfiV | cr | 

87rcr 2 Erfi 2 ,/R 




FIG. 3: Energy (in units of E a ) as a function of cr for various 
rj (N — 2k 3> 1). From top to bottom are curves for 77 = 9.9, 
9.96, 10, 10.0965, 10.3. 
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FIG. 4: (Color online) Nematic order parameter as a function 
of 77 for iV = 2fc(3> 1). Along the black lines are ground states; 
along the red (light) lines are metastable states. (See also the 
caption of FIG. 2.) 



in which Erfi[x] is the complex error function defined by 
Erf[za;]/z. In a series expansion for tr -C 1, the result is: 

zJcx ( 2 8 \ 2 

+ E s zJ ex <J A 

^315 14175 V 

( 8 32 \ 

-^ + 165375^J ff4 + ^) ^ 

The energy as a function of a at different 77 is plot- 
ted in FIG. 01 which is qualitatively the same as FIG. 
n for two particles per site. When rj < 9.96, the en- 
ergy as a function of a has only one (global) minimum, 
which corresponds to a spin singlet ground state. When 
77 > 9.96, in addition to the global minimum, there ap- 
pears a local minimum at a > 0. At 77 = rj c = 10.0965, 
these two minima become degenerate, signifying a phase 
transition. At rj > rj c , the solution at a = becomes a 
local minimum indicating a metastable spin singlet state, 
whereas the global minimum at a > corresponds to a 
nematic ground state. As 77 further increases, the solu- 
tion at a = becomes unstable and a local minimum 
occurs at a < 0, while the global minimum remains at 
a > 0. Following discussions on Eqs. 32,33 we inter- 
pret the a < solution as a metastable disk-like spin 
nematics. 

For the trial wave function in Eq. 44 the nematic order 



parameter can be calculated as: 

q = (cr\Q a /3W) 
(oo\Q a/ 3\oo) 



2 4cr 2 v / 7r|aiErfi v / |a| 

When Q is small, we obtain an expression of energy in 
terms of Q: 

(47) 

The jump in Q at the phase transition is equal to 0.323. 

The evolution of ground state wave functions and re- 
sults on quantum phase transitions are summarized in 
figure^l where the nematic order parameter is plotted as 
a function of 77. As stated before, these results are only 
valid in high dimensional lattices, where fluctuations in 
ordered states are small. For detailed calculations of fluc- 
tuations we refer to appendix iBl 



C. Large N limit: An odd number of particles per 

site 

At last, we also present results for an odd number of 
atoms per site. The main difference between this case and 
the case for an even number of particles per site is that at 
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zero hopping limit in the former case there is always an 
unpaired atom at each site. Consequently in the mean 
field approximation, we only find nematic Mott insulat- 
ing phases. As in the case for even numbers of particles 
per site, we expect this approximation to be valid in high 
dimensional lattices but fail in low dimensions, especially 
in one-dimensional lattices where long wave length fluc- 
tuations are substantial. Here we restrict ourselves to 
high dimensional lattices only. 

For large N a trial wave function which interpolates be- 
tween spin singlet states (dimerized) and nematic states 
can be introduced as: 

*odd({n fc }) = J] C{0,a) [0(n k ■ n )(n; • n ) + (n fe ■ n,)] 

(kt) p 

xexp[a((n fc -n ) 2 + (n r n n ) 2 )]. (48) 
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FIG. 5: The value of O as a function of 77 (N = 2k + 1 > 1). 



(fcZ)p denotes that the summation should be taken over 
parallely ordered pairs of nearest neighbors k and I cover- 
ing the lattice. C(0, a) is a normalization constant. The 
solution with O — 0, a — corresponds to a dimerized 
valence bond crystal state; and solutions with O 7^ 0, or 
o" ^ represent nematic states. 

It is straightforward, but tedious to compute the en- 
ergy of these states. Minimizing it with respect to various 
values of rj for d = 3 gives the results shown in FIG. 
and|^l No phase transitions arc found in the mean field 
approximation; and ground states break both rotational 
and translational symmetries^ 

At very small 77, the on-site Hilbert space is truncated 
into the one for a spin-one particle^ The reduced Hamil- 
tonian in the truncated space is a Bilinear-Biquadratic 
model for spin-1 lattices 




FIG. 6: The value of a as a function of r\ (N = 2k + 1 > 1). 



Hb.b. = J$>os0S fc • S, +sin#(S fc ■ S,) a ],Sj» = 2; (49) 
(kl) 

8 in general varies between — 3ir/4 and —ir/2. We there- 
for expect ground states at small rj limit should still ex- 
hibit nematic order (i.e. O ^ 0). 

It is worth emphasizing that conclusions about small 
rj limit arrived here are only valid in high dimensional 
bipartite lattices. In low dimensional lattices, states of 
correlated atoms in this limit were discussed recently and 
ground states could be rotationally invariant dimcrized- 
valence-bond crystals^ 



IV. CONCLUSIONS 

We have studied the microscopic wave functions of spin 
nematic and spin singlet Mott states. Both disk- like and 
rod-like spin nematic states were investigated. We also 
have analyzed quantum phase transitions between spin 
singlet Mott insulating states and nematic Mott insulat- 
ing states. We show that in the mean field approxima- 
tion, the phase transitions are weakly first order ones. 



Thus, we expect that fluctuations play a very important 
role in these transitions and the full theory on quantum 
phase transitions remains to be discovered. 

On the other hand, we have estimated fluctuations in 
different regimes of the parameter space. We found that 
fluctuations arc indeed small away from the critical point, 
at either small hopping or large hopping limit for an even 
number of particles per site. At the small hopping limit, 
fluctuations are proportional to 77, while at the large hop- 
ping limit they can be estimated to be proportional to 
(see Appendix B). 

For an odd number of particles per site, fluctuations 
are small only at large hopping limit and arc significant 
at small hopping limit. The later fact implies a large 
degeneracy of Mott states at zero hopping limit which 
was emphasized in the discussions on low dimensional 
Mott states. The physics in this limit remains to be fully 
understood. 

In the context of antifcrromagnets, spin nematic states 
have also been proposed plSLSSL^i Collective excitations in 
atomic nematic states should be similar to those studied 
in previous works; we present some brief discussions on 
this subject in Appendix B and refer toi2*2£L2I for details. 
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APPENDIX B: SPIN FLUCTUATIONS IN MOTT 
STATES 

Nonlinear dynamics and spin fluctuations in conden- 
sates of spin-one bosons were discussed in a previous 
work— Here we carry out a similar discussion for Mott 
states. Following the Hamiltonian 

H = E s J2^l~ Jcx]T(n fc -n ; ) 2 (Bl) 

k (kl) 

we derive the following equation of motion for the direc- 
tor rife: 



APPENDIX A: AN ALTERNATIVE 
DESCRIPTION 

Alternatively, one can also carry out the calculations 
in section III, using the following operator: 



?2,aa' = S Q S Q — — S aa i S 7 S 7 
_ _,a07,a'0y,/.t,/.t 



(Al) 



—— = 2E s S k x n fe 
at 



(B2) 



Fluctuations when 77 is small 

For rj = zJ ex / E s = and an even number of particles 
per site, the ground state is the product state: 



tf„=o = n y oo(n fc ) 



(B3) 



When <)) < 1, the ground state wave function can 
also be obtained by a perturbation theory; the leading 
term is 



Defining the more conventional order parameter— 

/=. \Q2,aa') 



(Q2,aa' 



(A2) 



ref 



we obtain the following results for the trial wave function 
in Eq. 28: 



Q-2 

E 



6E t 



sin 2 6 



lzJ ex I 2^21 




28 - ~ 
—zJ cx Q 2 2 + 0(Q 



5/2, 



which lead to the same conclusions as in section III. How- 
ever, in terms of the order parameter defined in Eq. A2, 
the rod-like and disk-like structures shown in FIG. 2 and 
FIG. 4 are less obvious. 

In the case of a large number of particles per site, the 
order parameter introduced here has the same expecta- 
tion value as the operator in Eq. 40. 



* ( o ] = E 



P (0) p(0) 



(B4) 



In our case = Y/ m with I even, and E[ = l(l + l)E s 
A direct calculation yields 



(i) _ JO 
~ 45 



;E E Y 2m (n t )Y 2 ^ m (n,) JJ r 00 (n fe ). 



(ij) m=-2 



(B5) 

Taking into account (ioolQa/sl^oo) = 0, wc find desired 
results in this limit, 



(Qk,ap) — 0. 



(B6) 



To characterize fluctuations, we study the following 
correlation function (Qk.aaQk' .aa}- Calculations of this 
correlation function in the state given in Eq. IB5I yield 

(Qk,aaQk',aa) = ~rpS{kk' , (kl))x 

45 



'E ( (^m|Qaa|5^o)(^,-m|Qaa|^Oo) + h.i 



S(kk', (kl)) is unity if k! and k sites are two neighboring 
sites as (kl) and otherwise is zero. The last expression 
can be calculated explicitly, 



(V2 in \Qaa\Yoo)(Y2,- m \Qaa\Yoo) + 



1 = 45 (B7) 
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Clearly at small 77, fluctuations are small. 
Fluctuations when 77 is large 

Again we consider the case for an even number of par- 
ticles per site. In the limit of rj — > 00, all directors n^ 
point in the direction of e 2 . For a finite but large rj we 
introduce 



i - r 2 

-L L 1, 



kx 



c 



(B8) 



where C^q,, a = .t, y are much less than unity. 

Following discussions in section IIIB, we obtain the 
following commutators. 

[Sky,Ck' x ] ~ i8k,k', [Skx,Ck> y ] ~ —iSk,k' (B9) 
which define two sets of harmonic oscillators. Introducing 



To obtain cigcnmodes, we perform a Fourier transfor- 
mation (setting the lattice spacing to be unity), 



rr \ yr Z7rfc-q s~1 _ \ ^ ( < Z7r/c-q 



(B13) 

where Vt is the total number of lattice sites. This leads 
to the following Hamiltonian 



(B14) 



q,a 



Following a standard calculation, fluctuations in this 
limit are: 



the effective Hamiltonian becomes 



(BIO) 



a=x.y k (kl) 



And 



[n^c, Cfc/^] = iSk y k'S a 0- 



(B12) 



- E 



\^?^ | q ^ c sin|q|7r/2 
(B15) 

The momentum cut-off q c in general depends on the short 
distance behavior of our model and for simplicity we set 
it as one. In high dimensional lattices, the sum in Eq. 
IB 151 is convergent; and we see the fluctuations are also 
small at the large rj limit. 
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